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Abstract 

Directional  wave  energy  and  wave  spectra 
play  an  important  role  in  the  physical  processes 
associated  with  the  ocean  environment. 
Determining  these  directional  wave  characteristics 
is  a  fairly  arduous  task.  An  Autonomous 
Undem’ater  Vehicle  (AUV)  provides  a  suitable 
platform  for  making  in-situ  measurements  in  the 
ocean  environment.  Small  AUVs  can  typically  cover 
a  3-5  mile  square  region  at  a  speed  of  3-4  knots, 
allowing  for  the  survey  of  the  water  column  over  a 
substantial  range.  Further,  an  AUV  is  much  more 
versatile  than  the  traditional  bottom-mounted 
ADCPs,  suspended  current  meters  or  drifting  buoys 
in  so  much  as  the  sensor  (the  AUV)  is  mobile  and 
can  be  programmed  to  collect  data  in  any  desired 
area  vice  having  to  locate  the  sensors,  recover 
them,  transit  and  re-deploy. 

In  this  paper  we  examine  the  potential  of 
an  AUV  to  be  used  as  a  platform  for  making 
directional  wave  spectrum  measurements,  with 
particular  reference  to  the  Naval  Postgraduate 
School’s  PHOENIX  AUV.  We  show  that  by  using 
relative  velocity  measurements  from  a  SONTEK 
Acoustic  Doppler  Velocimeter  (ADV),  ground- 
referenced  velocity  measurements  from  a  RDI 
Navigator  Doppler  Velocity  Log  (DVL),  vehicle 
motion  data  from  a  Systron-Donner  Motion  Pak 
and  employing  a  Maximum  Entropy  Method 
(MEM)  the  wave  directional  spectra  can  be  found. 
The  methodology  used  to  determine  the  spectra,  the 
corrections  that  are  required  to  account  for  the 
Doppler  shift  due  to  the  moving  vehicle  and  the 
results  obtained  from  data  collected  during 
AUVFEST  ‘98  and  experimental  missions  in 
Monterey  Bay  will  be  discussed  and  analyzed. 

Introduction 

Measurements  of  directional  wave  spectra  have 
interests  that  range  from  fundamental  physics,  to 
accurate  forecasting,  to  determination  of  tactical 
information  for  military  operations.  There  are 
many  techniques  for  collecting  spectral 


information,  each  of  which  has  advantages  and 
disadvantages.  These  include  deployable  or  in-situ 
instruments  such  as  directional  wave  buoys 
[Grosskopf  1983],  as  well  as  arrays  of  current 
meters  and  pressure  sensors  [Allender  1989]. 
Remote  sensing  techniques  involving  microwave 
radar  systems  [Tyler  1974],  aircraft  [McLeish 
1980],  or  satellites  [Monaldo  1984]  are  also 
commonly  used.  The  issue  associated  with  these 
measuring  approaches  is  that  they  are  too  expensive 
for  many  users  and  are  not  suitable  for  rapid 
deployment  without  significant  pre-planning. 

Currently  many  operational  AUVs  carry 
sensor  suites,  which  allow  the  vehicles  to  record 
data  that  many  be  used  to  obtain  directional  spectra 
estimates.  This  paper  will  outline  the  underlying 
principles  used  in  identification  of  wave  directions 
from  standard  wave  following  buoys.  It  will 
present  the  mathematical  formulas  used  in 
determining  the  wave  direction  as  a  function  of 
frequency.  Extension  of  these  algorithms  to 
subsurface  velocity  sensors  will  be  made,  where, 
through  the  use  of  the  Doppler  equation,  a  moving 
AUV  can  be  used  to  determine  wave  directions. 
Lastly,  it  will  be  shown  how  a  control  command 
can  be  obtained  from  the  frequency  dependent 
wave  direction  estimates. 

The  information  in  this  paper  is  not  new, 
only  the  application  to  which  this  method  is 
applied.  For  more  detailed  descriptions  of  the 
mathematical  formulations  presented  in  this  paper, 
the  reader  is  referred  to  [O'Reilly  1996]  and  the 
references  therein. 

Wave  Spectra  and  Directional  Estimates 

The  elevation  of  the  sea  surface  can 

be  described  as  the  superposition  of  an  infinite 
number  of  sinusoids  of  the  form: 

iiW=X>„  c°s(m-(v+cu=2X  ■  (!) 

n= 1  n= 1 

Wave  components  with  different  frequencies  are 
usually  assumed  to  be  statistically  independent 
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because  they  are  generated  by  random  wind  forces 
at  different  locations.  From  the  central  limit 
theorem  it  follows  that  the  probability  distribution 
of  the  surface  elevation,  r\(t),  is  approximately 
Gaussian,  consistent  with  many  observations, 
[Soong  1993], 

The  procedure  presently  employed  by 
many  of  the  operational  data  buoys  in  based  on 
Fourier  analysis.  In  Fourier  analysis  it  is 
convenient  to  work  with  complex  exponentials 
rather  than  sine  and  cosine  functions,  there  fore 
using  the  relation 

cos(0)f  +4) )  =  exp(i(m ^  +  CXp(-,tof  )}  (2) 

the  expression  for  the  surface  elevation  can  be 
written  as 

!](?)  =  EA»  exp(icor) ,  (3) 

0 

where 

exP('40, )  (4) 

and  the  summation  is  over  both  positive  and 
negative  frequencies. 

The  energy  spectrum  E( to),  is  defined  as 


where  (  indicates  an  average  over  many  data 


records  and  Aco  is  the  spacing  of  the  frequency 
bands.  The  spectrum  is  closely  related  to  the 
energy  of  the  waves,  and  represents  the  distribution 
of  wave  energy  as  a  function  of  frequency. 

To  describe  the  spatial  and  temporal 
characteristics  of  the  sea  surface  linear 
superposition  of  wave  components  is  used.  In 
exponential  terms  this  can  be  represented  as 
q(v,  y,t)  =  ILAo  e  exp(/(£(xcos0  +  ysinO)  —  wt)) 

CO  0 


(6) 

where  x,  y  are  the  horizontal  spatial  coordinates, 
and  to  and  k  obey  the  dispersion  relation.  The 
frequency  directional  wave  spectrum  is  defined  as 


£(©)  = 


^o),e 


Aco  A0 


(7) 


and  describes  the  distribution  of  energy  as  a 
function  of  frequency  and  direction. 


Directional  Estimates  from  Wave  Buoys 


The  most  commonly  used  instrument  for 
measuring  waves  in  deep  water  is  the  "heave,  pitch 
and  roll  buoy"  that  measures  the  surface  height  and 
slope  in  two  orthogonal  directions.  The  newer 


Datawell®  Directional  Waverider  measures 
3-component  accelerations  of  the  buoy,  which  are 
integrated  to  yield  the  horizontal  and  vertical 
displacements  of  the  buoy.  The  hull  and  mooring, 
of  this  buoy,  were  designed  to  give  the  buoy  good 
wave  following  characteristics,  thereby  allowing 
the  buoy  displacements  to  approximate  the 
displacements  of  an  actual  water  particle  at  the  sea 
surface. 

For  a  wave  train  propagating  in  the 
positive  x-direction,  the  fluid  particle  motion  is 
given  by 

gk  cosh(k(z  +  H))  , 

u  =  a- - cos(fcr-cof) .  (8) 

co  cosh  (kH) 

However,  for  the  more  general  case  of  a  wave  train 
propagating  at  some  angle  relative  to  the  x-axis,  it 
can  be  shown  that  the  flow  field  is  given  by 
m(x,  y,t )  =  acocosO  exp(LZ)cos(£(xcosO  +  y  sin0)  -  cot) 

v(x, y,f)  =  aco  sinO  exp(AZ)cos(£(xcos0  +  vsinO)— cot) 
w(x,  yj)  =  aco  exp(fcZ)sin(A:(xcos0  +  y  sinO)  —  (Of) 

(9) 

Let  the  average  position  of  the  wave  buoy  be  given 
by  x=y=z= 0.  For  small  amplitude  waves,  the 
expected  buoy  displacements  are  small  compared  to 
the  surface  wavelength,  therefore  the  buoy  motion 
can  be  approximated  by  the  fluid  velocity  at 
x=y=z= 0. 

For  a  full  spectrum  of  waves,  the  buoy 
displacements  can  be  expressed  using  complex 
notation  as 

x ( t )  =  EE_*A,>.e  cos6  exp (iwt) 

0  0 

Y(t)  =  EE_/Ao.e  sin0  exp (iwt)  (10) 

0  0 

m  =  EE-4,,,0  exp  (iwt) 

0  0 

where  the  -i  is  due  to  the  90  phase  difference 
between  the  vertical  and  horizontal  displacements. 

The  expressions  in  (10)  can  be  written  using 
Fourier  transforms  as 

X  ( t )  =  X  (co)exp(iwf) 

0 

Y(t)  =  ^Y((X>)exp(iwt)  ,  (11) 

0 

Z(t)  =  ^Z(co)exp(wf) 

0 

where  the  Fourier  transforms  are  given  by 

Z  (co )  =  E  —iA<ng  cos0 
0 

Y (co)  =  E  —  sin0  ■  (12) 

0 

Z((0)  =  EA, 9 

e 
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To  derive  the  relationships  between  the 
measured  time  series  and  the  unknown 


frequency-directional  wave  spectrum  the  cross 
spectrum  must  be  considered.  In  general,  the  cross 
spectrum  between  two  time  series  X(t)  and  Y(t )  with 
Fourier  transforms  AY  op  and  Y(  co)  is  defined  as 


CXY(  co) 


(X  (00)7*  (to)) 
Aco 


(13) 

where  indicates  the  complex  conjugate,  [Soong 
1993],  Substitution  of  (12)  into  (13)  yields 

CXY  (co)  =  J^cosG  sinGEYco.G ) 

e 


(14) 

where  it  is  assumed  that  the  wave  components 
propagating  in  different  directions  are  statistically 
independent.  The  cross  spectrum  CXy  can  be 
expressed  in  continuous  form  as 

271 

Cxy( co)=  JcosG  sinG£(co,6)c/6  . 

o 

(15) 

Cross-spectra  of  the  other  time  series  pairs  can  be 
obtained  in  a  similar  manner.  The  full  set  of 
relations  for  buoy  cross-spectra  can  be  found  in 
[Dean  1984], 

It  is  convenient  to  define  a  normalized 
directional  distribution  of  energy  at  a  frequency  co 
as 


S(G;co)  = 


E(co,6) 
E(  co) 


(16) 


with  unit  integral 


JS(G;coV/G  =  — 


j£(co,G)dG 


E(  co) 


E(a) 
E(  co) 


=  1. 


(17) 

With  this  definition,  (15)  and  the  other  referenced 
spectral  relations  can  be  combined  and  expressed  in 
terms  of  5(0  ;  CO )  .  Dropping  the  frequency 
dependence  these  relations  can  be  expresses  as 
Im  (Cxz) 


(((-  XX  +  Cyy  )czz  ) 


=  JcosGS(G  WG  =  flj 


Im(CK) 


((CXX  +  Cyy  )CZZ  ) 


r  -C 

^XX  ^YY 

c  +r 

^yy 


271 

JsinG5(G)c/G  =bt 

0 


271 

Jcos2G5(G)c/G  =  a2 
o 


2R  e(CXY) 
C  +C 

^XX  T  ^YY 


2k 

=  Jsin  2G.SYG  )c/0  =  b2 

0 


(18) 

These  four  relations  between  the  cross-spectra  of 
the  buoy  measurements  and  the  directional 
distribution  of  wave  energy,  derived  by  [Long 


1980],  form  the  basis  for  most  of  the  buoy  analysis 
techniques  currently  used. 

Extension  to  Subsurface  Sensors 


As  discussed  in  the  previous  section,  the 
motion  of  a  wave  buoy  is  directly  related  to  the 
fluid  velocity,  therefore,  the  cross-spectra  of  a  tri- 
directional  current  meter  yields  the  same  low- 
resolution  directional  wave  information  obtained 
from  buoy  measurements.  Substituting  the 
normalized  spectra  of  the  vertical  (Z)  velocity 
component  w,  and  the  horizontal  (X,  Y)  velocity 
components  u  and  v  into  (18),  the  lowest  four 
Fourier  moments  of  the  directional  distribution  of 
wave  energy  can  be  obtained  and  are  given  by 
Im(C  (co)) 


a,(co)  = 


bt(  «)  = 


a2(co)  = 


M  ®)  = 


[Cw({o)[c„„(co)  +  Cvv(o))]'/ 
Im(COT(co)) 

[cw(co)[c„„(co)  +  Cw(®)F 

cinim-cvvm 

C„„(co)  +  Cw(co) ' 

2Re(C„v(co)) 


(19) 


(20) 


(21) 


(22) 


C„„(co)  +  C„,(co) 
where  C((0)  is  the  spectral  matrix  of  the  velocity 
components  u,  v,  w.  Since  the  wave  direction,  G,  is 
referenced  to  the  navigation  frame  (N-E-D),  vehicle 
borne  sensor  measurements  must  be  transformed 
prior  to  use.  It  is  interesting  to  note  that  the 
estimates  of  these  directional  moments  are 
insensitive  to  errors,  so  long  as  the  errors  are  the 
same  on  all  measurement  axes  of  the  sensors,  which 
is  typical  with  oceanographic  sensors  installed  on 
AUVs. 

The  objective  of  the  data  analysis  is  to 
infer  the  directional  distribution  5(G),  from  the  four 
measured  moments  ah  bh  a2,  and  b2.  The  most 
widely  used  techniques  are  described  below. 


The  Cosine  Power  Distribution 


[Longuet-Higgins  1963]  introduced  a 
simple  cosine-power  distribution, 

5(6)  =  Acos2^9  j  (23) 

with  G, the  mean  propagation  direction,  ,v  a 
parameter  that  controls  the  width  of  the  distribution 
and  A ,  a  normalization  coefficient.  The  parameters 
Qmean  and  .v  are  determined  by  fitting  (23)  to  the 
relations  given  in  Equations  (19)-(22). 

The  main  drawback  to  this  simple  method  is  that 
(23),  with  only  two  free  parameters  can  describe 


3 


only  unimodal  distributions,  and  thus  fails  in 
situations  with  a  bimodal  sea  state  (e.g.,  multi¬ 
directional  seas  during  a  wind  veering  event  or 
swell  arriving  from  two  different  sources). 


The  Maximum  Entropy  Method 


[Lygre  and  Krogstad  1986]  introduced  the 
maximum  entropy  or  MEM  estimate  of  S( 0). 
Unlike  Equation  (23),  this  approach  can  describe 
both  unimodal  and  bimodal  distributions  and 
exactly  fits  the  relations  given  in  (19)-(22).  This 
directional  spectrum  is  given  by 


i  i— <>ic-r  —<>2^2 

2n  (l-t^e  ,e  — 2,0  )~ 


(24) 


with 


q  =  aI  +  ibt 


c2  =  a0  +  ib2 


<l>i  = 


Ci  c,c. 


1-  c, 


(25) 


<l>2  =c2-ct( j), 

Still,  the  directional  distribution  is  poorly 
constrained  by  only  four  moments  and  the  estimates 
require  careful  interpretation,  [Krogstad  1991], 


Mean  Direction  and  Directional  Spread 


An  alternative  approach  that  avoids  the 
pitfalls  of  5(0)  estimation,  is  to  describe  the 
directionality  of  waves  by  a  few  simple  parameters. 
For  narrow  5(0 ),  a  mean  propagation  direction  0m 
and  a  root-mean-square  measure  of  the  directional 


spreading  energy  O0  can  be  defined  in  terms  of 
the  first-order  and  second-order  Fourier  moments 
a,,  bx,  a 2  and  b2  [Kuik  1988],  given  by, 

(u  \ 


9,„  =  tan 


h 


(26) 


a2  =  2[l  -  [a,  cos(0J  +  bx  sin(0,„ )]  (27) 

1 


e,„=-tan" 
1 


fb^ 


(28) 


oe2  =  — [l  —  [«2  cos(20m)  +  /?2  sin(20m )].  (29) 


Again,  this  method  fails  to  identify  bimodal 
distributions  but  it  is  useful  to  determine  a  base 
direction  so  that  a  control  command  could  be 
determined.  More  on  this  approach  will  be 
discussed  later  in  this  chapter. 


Correction  For  A  Moving  Platform  (The 
Doppler  Equation) 

The  equations  for  the  wave  directionality 
estimations  presented  in  the  previous  section  are 
valid  for  a  non-moving  sensor.  However,  when 
information  is  obtained  from  an  AUV,  corrections 
must  be  made  to  account  for  the  vehicle  motion. 
The  wave  frequency,  which  the  vehicle  encounters 
while  moving  through  the  wave  field,  has  been 
shifted.  This  shift  can  be  determined  by  using  the 
well-known  Doppler  equation, 

2 

co^  =co - U  cos  p  ,  (30) 

5 

where  the  spectrum  that  the  vehicle  encounters  is  a 
function  of  the  vehicle's  forward  speed  U  and  its 
heading  angle  relative  to  the  propagation  direction 
of  the  sea  waves,  [j.  Using  the  techniques  outlined 
in  the  previous  section  will  give  the  wave 
directional  distribution  as  a  function  of  vehicle 
encounter  frequency.  If  these  estimations  are  used 
in  conjunction  with  the  Doppler  equation  in  a 
recursive  manner,  the  estimation  of  5(0)  can  be 
found  as  a  function  of  true  frequency. 

The  method  by  which  this  is  performed  is 
outlined  below; 

1.  Determine  the  three  components  of  the  fluid 
velocity  in  vehicle  body  fixed  coordinates. 

2.  Transform  the  fluid  velocities  into  the  global 
navigation  frame. 

3.  Compute  the  auto-  and  cross-spectra  of  the 
fluid  velocity  components. 

4.  Determine  the  Fourier  moments  using  (19)- 

(22). 

5.  Convert  the  Fourier  moments  into  Krogstad 
notation  and  use  the  MEM  to  determine  the 
directional  distribution  5(0). 

6.  Using  the  vehicles  mean  velocity,  and  the 
mean  direction  obtained  from  (26)  or  (28), 
apply  the  Doppler  equation  to  determine  the 
frequency  shift. 

7.  Return  to  3  and  complete  the  process  until 
frequencies  converge. 

The  corrections  to  the  estimation  of  5(0)  using  the 
Doppler  equation  are  quite  small  for  slow  moving 
vehicles.  Considering,  for  example,  the  NPS 
Phoenix  AUV  which  has  a  maximum  forward 
velocity  of  1.5  m/s,  the  error  associated  with  not 
using  the  Doppler  equation  are  approximately  ±  1 
sec  in  identification  of  wave  periods  between  4  and 
40  seconds.  Similarly,  the  error  in  direction 
estimation  is  approximately  5-7  degrees.  When  an 
AUV  goes  into  a  station-keeping  mode,  where 
vehicle  velocities  are  significantly  reduced,  the 


4 


modifications  required  due  to  the  Doppler  shift  are 
negligible. 

Experimental  Results  Of  Wave  Directional 
Estimation  Using  The  NPS  Phoenix  AUV 

During  November  1998,  the  NPS  Center 
for  AUV  Research,  under  the  direction  of  Professor 
Anthony  Healey,  participated  in  the  ONR 

sponsored  AUVFEST  '98.  This  AUV  technology 
demonstrator  was  held  in  the  Gulf  of  Mexico,  south 
of  Gulfport,  MS.  A  complete  description  of  the 
exercise  can  be  found  in  [Bunce  1998], 

The  Phoenix  AUV  was  used  during  this  exercise  as 
a  mobile  sensor  to  gather  oceanographic  data. 
Using  the  concepts  presented  in  the  previous 
section,  the  vehicle  conducted  27  short-term 

sampling  missions.  The  products  that  were 
obtained  included  directional  energy  plots, 

directional  spectrum  plots  and  mean  current 

estimations. 

The  key  to  the  ability  of  Phoenix  to  obtain 
data  capable  of  producing  this  information  is  the 
combined  ADV/RDI/MotionPak/Directional  Gyro 
sensor  suite.  The  physical  layout  of  the  NPS 
Phoenix  AUV  is  shown  in  Figure  1.  Detailed 
description  of  the  vehicle  can  be  found  in  [Marco 
1996]  and  [Brutzman  1997].  In  addition  an  online 
description  can  be  found  on  the  NPS  center  for 
AUV  Research  web  site  at 
http://www.cs.nps.navy.mil/research/auv/auvstats. 
With  these  sensors,  accurate  three-dimensional 
fluid  velocities,  expressed  in  global  quantities,  were 
capable  of  being  obtained  in  post-processing.  Since 
the  RDI/ADV  sensors  are  collocated,  little  vehicle 
induced  motion  remains  after  processing  the  data. 

The  data  obtained  validated  the  concept  of 
obtaining  tactical  oceanographic  data  from  an 
underwater  vehicle.  During  the  collection  of  the 
data,  remnants  of  Hurricane  Mitch  were  still 
present  in  the  Gulf,  providing  an  excellent  source  of 
ground  swell.  In  addition,  there  was  a  significant 
wind  generated  wave  component  in  a  different 
direction  than  the  swell  component,  resulting  in  a 
multi-modal  spectrum.  Figure  2. 

The  results  presented  in  Figures  3-5 
provide  a  sample  of  the  oceanographic  data 


obtained  during  this  offshore  exercise.  As  seen  in 
Figure  3,  the  bi-modal  properties  of  the  seaway  are 
captured,  as  well  as  an  estimate  of  the  significant 
wave  height.  The  ability  to  estimate  the  dual 
directions  is  due  to  the  use  of  the  MEM  algorithm 
employed.  Figure  4  presents  the  associated 
spectrum  plots  for  this  energy  estimate.  The  ability 
of  an  AUV  to  estimate  currents  is  shown  in  Figure 
5.  Using  a  triangular,  time  based  run;  the  current 
can  be  determined  using  set  and  drift  calculations 
from  the  error  in  final  position  as  well  as  the 
heading  error  on  each  leg.  This  information  can  be 
feed  directly  into  the  vehicle’s  navigation  process 
to  account  for  the  offset  due  to  current  thereby 
increasing  navigation  accuracy.  Short-term 
averages  from  each  of  the  three  legs  are  in  general 
agreement  with  the  overall  average  determined 
from  the  navigational  drift. 

During  experimental  validation  of  a  new 
Disturbance  Compensation  Controller  [Riedel 
1999],  wave  direction  estimates  were  obtained  in 
Monterey  Bay,  see  Figure  6.  Again  the  Phoenix 
was  able  to  record  data  which  can  be  used  to 
determine  directionality  spectra.  Figure  7. 

Conclusions 

This  paper  has  presented  the  techniques 
currently  employed  for  the  determination  of  wave 
directional  estimations  from  standard  wave 
following  buoys.  It  has  extended  this  analysis  for 
use  in  determining  directional  estimates  from  data 
gathered  by  an  Autonomous  Underwater  Vehicle. 
Using  this  data  gathered,  an  approach  was 
presented  which  will  allow  the  deployed  vehicle  to 
obtain  information  about  the  directionality  of  its 
working  environment  thereby  allowing  the  vehicle 
to  have  information  available  to  make  decisions 
regarding  mission  execution.  Actual  experimental 
results  using  the  NPS  Phoenix  AUV  have  been 
presented.  It  was  shown  that  tactical  oceanographic 
data  is  obtainable  from  a  moving  AUV.  The 
vehicle  in  this  manner  becomes  an  intelligent, 
mobile  off-board  sensor,  thus  presenting  the 
argument  for  AUV  deployment  with  operational 
fleet  units. 
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Figure  1  Physical  Layout  of  Phoenix  AUV 


Figure  2  Phoenix  Being  Launched  from  the  R/V  Gyre 
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Phoenix  Survey  Data 
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Figure  3  Sample  Direction  Energy  Plot  From  Phoenix  Data,  Nov.  4,  1998  (Run*  9),  Gulf  of  Mexico, 

AUVFEST  ‘98 
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Figure  4  Sample  Direction  Spectrum  Plots  From  Phoenix  Data,  Nov.  4,  1998  (Run  9),  Gulf  of  Mexico, 

AUVFEST  ‘98 
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Figure  5  Short-term  Current  Estimation  from  Phoenix,  November  8,  1999,  (Run#  2),  Gulf  of  Mexico, 

AUVFEST  ‘98 
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Figure  6  Sample  Direction  Energy  Plot  From  Phoenix  Data,  June  15,  1999  (Run#  4),  Monterey  Bay 
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Figure  7  Sample  Direction  Spectrum  Plots  From  Phoenix  Data,  June  15,1999  (Run*  4),  Monterey  Bay 
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